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Abstract 

Based on the permutation group formalism, we present a discrete symmetry algebra in QCD. The 
discrete algebra is hidden symmetry in QCD, which is manifest only on a space-manifold with non- 
trivial topology. Quark confinement in the presence of the dynamical quarks is discussed in terms 
of the discrete symmetry algebra. It is shown that the quark deconfinement phase has the ground 
state degeneracy depending on the topology of the space, which gives a gauge-invariant distinction 
between the confinement and deconfinement phases. We also point out that new quantum numbers 
relating to the fractional quantum Hall effect exist in the deconfinement phase. 

PACS numbers: 03.65 Fd, 12.38 -t, 12.38 Aw 
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The purpose of this paper is to present an argument for classification of phases in QCD. 



Classification of phases in QCD is an old but unsolved problem in quantum field theory. 
(For recent reviews, see [l], As is well-known, behaviors of the Wilson loop [2I (or 
Polyakov line {4]) and the 't Hooft loop jl], 6] are useful to classify the quark confinement 
and deconfinement phases in the pure Yang Mills theories, but once the dynamical quarks 
are included, they are no longer sufficient to distinguish them. Nevertheless, it will be shown 
below that there exist quantum numbers which distinguishes these two phases. In particular, 
they are useful to study the quark confinement in the numerical studies. 

This work is motivated by recent developments on understanding of quantum phases. 
Many quantum phases and phase transitions can be described by the spontaneous symmetry 
breaking and local order parameters, however in recent years it has become increasingly clear 
that in a wide class of strongly correlated many-body systems, a phase transition driven by 
a nonthermal parameter may occur at zero temperature which can not be understood by 
any local order parameter. The characteristic signature of the novel phase is a finite ground- 
state degeneracy depending on the topology of the space, and the underling order of the 
novel phase is dubbed as topological order j^K The Laughlin state for the fractional Hall 



effect is known to have a topological order 8 



In the present, many systems including 



josonic ones and those at zero magnetic field are identified as possessing topological orders 
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Recently, we have argued that the topological degeneracy in a topological order is due 
to the emergence of a discrete symmetry [18], which contains three fractional parameters: 
quasiparticle charge, anyon statistics, and the fractional quantum Hall conductivity. In par- 
ticular, it is notable that the emergence of collective excitations having fractional quantum 
numbers with respect to the constitute particles in the Hamiltonian is closely related to the 



existence of the topological order 18j, ll9[. Such a charge fractionalization has an interesting 



analogy of the quark deconfinement, which also gives fractional charged excitations. In spite 
of the essential difference that the quarks are elementary particles, not collective excitations, 
this motivates us to study the quark confinement in the notion of the topological order. 

In the following, we will show that the quark deconfinement phase in QCD indeed has a 
topological order. Generalizing the argument in 18j, we will construct a discrete symmetry 
algebra in QCD, which we dub topological discrete algebra. The existence of the center of 
the gauge group is crucial for the construction. By the use of the topological discrete algebra, 
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it will be shown that the quark deconfinement phase in QCD has a ground state degeneracy 
depending on the topology of the system. The topological degeneracy in the thermodynamic 
limit is a gauge-invariant quantum number that distinguishes the deconfinement phase from 
the confinement one. 

For definiteness we will consider the lattice QCD. The generalization to the continuum 
QCD is straightforward. The action of the link variable U n>ll G SU (3) is given by 

p ^ 

with the plaquette variable Up = C/ n ,/i^n+A^^n+i> uPn,vi an d that of the quarks -0™ is 

Sf = -^EE (feMi/i - €^ui^Cj - E m rfM, (2) 

/ n,fj, f,n 

where n = (jii, n 2 , n%, n^) denotes the site on the lattice, fx a unit vector in the direction, 
and / the indices of the flavors of the quarks. To remove the doublers, we also add the 
Wilson term S\v- The partition function Z is given by 

Z = J VUVtpVije- {SG+SF+Sw) . (3) 

If necessary, other terms preserving the electric charge (improved actions for gauge field, the 
chemical potential terms for quarks, and so on) can be included, which does not affect the 
arguments in the following. 

The topological discrete algebra we will consider is defined only if the topology of the 
space-manifold is non-trivial, so let us consider the system on a three dimensional torus 
T 3 . (The space-time is four dimensional.) On the torus T 3 , the site at n — (n 1; n 2 , n 3 ) 
is identified with one at n + aN a , where a is a unit vector in the n a direction, and iV a is 
an integer (a = 1,2,3). The torus size is NiN 2 N 3 . In practice, the torus is realized by 
imposing the periodic boundary conditions in all the spatial directions. The torus admits 
three independent spatial non-contractable loops C a (a — 1, 2, 3) along the n a direction. In 
addition, the torus has three "holes" h a (a = 1, 2, 3) outside T 3 . The hole h a is encircled 
by the non-contractable loop C a . In other words, each S 1 in T 3 = S 1 x S 1 x S 1 has a 
noncontractable loop and a hole. See FigJH 

To define the topological discrete algebra, we introduce an external U(l) electro-magnetic 
gauge field e ldn -^. Then consider an adiabatic insertion of the flux $ a through "hole" h a 
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FIG. 1: Non-contractable loop C a on S . 

encircled by the non-contractable loop C a . This process induces the external U(l) electro- 
magnetic gauge field e ldn ^ with 9 n ^ = 5^ a Q a /N a . Note that the £7(1) field strength on T 3 
remains to be zero after the flux insertion. 

Let us now consider the partition function 2>($> a ) with the inserted flux $ . Since the 
upper and down quarks have 2/3 and —1/3 electric charges, respectively, then by the unitary 
transformation, 

ip n -> e~ i2na ' s ' a/3Na ip n , for upper quarks, 

ip Tl -> e ina ^ a,3Na ilj ni for down quarks, (4) 

the induced U(l) electromagnetic gauge field is eliminated in the action except on the links 
connecting the sites n with n a = 1 and n a = N a . After the transformation the kinetic 
terms of the quarks on these boundaries acquire the U(l) phase as 

t l2 ^ atZ ^nlaV n ^n+a, for upper quarks, 

e~^ a/3 'ip n j a U nia ijj n+& , for down quarks, (n a = N a ). (5) 

The same U(l) phase remains in the kinetic term in the Wilson term. If $ a is a unit flux 27r, 
these U(l) phases coincide with an element of the center of SU(3), e -27 ™/ 3 . So performing 
the unitary transformation on link variables U nA on the boundaries, U U;a e 2nl ^ 3 U nA , (n a = 
N a ), one can show Z(2ir) = Z(0). Therefore, the spectrum of the system is invariant under 



the adiabatic flux insertion by 2ir. The unit flux 2ir is smaller than that in [18( where the 
unit flux is 2n/e with e the minimal charge of the constituent particles. In QCD, e = 1/3. 
The adiabatic insertion of the unit flux 2ir is represented by a unitary operator U a . 

If quarks are deconfined, the physical states are classified by the representation of the 
permutation group for quarks. For iV quarks on T 3 , the permutation group consists of <7j 



(z = 1, • • • , N— 1), which exchanges the ith and (i + l)th quarks clockwise without enclosing 
any other quark, and rf (a = 1, 2, 3, i — 1, ■ ■ ■ , N), which represents moving the ith quark 
along the non-contractable loop C a in the n a direction. The permutation group is given by 

a 2 k = l, l<k<N-l, 
{a k a k+1 ) 3 = 1, 1 < k < N - 2, 
o k oi = ai<Jk, 1 < k < N — 3, \l — k\ > 2, 
ri+i = vtfvi, 1 < % < N - 1, a = 1, 2, 3, 
^ = ^r", 2<j<N, a = 1, 2, 3, 
r " r i= r i r ^ U = l, a, 6 = 1,2, 3. (6) 

The generators of the permutation group have non-trivial commutation relations with U a 
(a = 1, 2, 3). If we apply rf after the flux insertion U a , the gauge field e lBn ^ will give rise to 
an Aharanov-Bohm phase e~ l2n ^ 3 . (Both the upper and down quarks acquire the same U(l) 
phase.) Therefore, we obtain 

jfU a = e- 2m / 3 U a r?, (7) 

where a = 1, 2, 3 and i — 1, • • • , N. On the other hand, because r\ (b ^ a) and cr^ do not 
encircle the inserted flux by U a , they commute with U a , 

rtU a = U a rl a,U a = U a a t . (8) 

with a 7^ b. 

Using the commutation relations (JTj) and ([8]), one can verify that UaUbU^U^ 1 (a ^ b) 
commutes with all the permutation group generators. So by Schur's lemma, for any irre- 
ducible representation of the permutation group, UJJ b U~ x U^ x is a (unimodular) c-number. 
Namely, 

U a U b = e 2mX ^U b U a , A 0i & = -A M . (9) 

In addition, U 3 is also a unimodular integer because it also commutes with all the permuta- 
tion group generators. Therefore U 3 U b = U b U 3 . Comparing this with we find that X a ^ b 
should be a rational number 

K b = k -f, (io) 
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where k a ^ is an integer and co-prime to 3. These constants X a>b {a,b = 1,2,3) are new 
quantum numbers in the deconfmement phase. 

If the system is time-reversal invariant, these new quantum numbers are shown to be zero: 
For the time-reversal invariant system, the time-reversal transformation of U a is given by 
TU a T~ l = CfJJ\ with a unimodular constant c a . Applying the time-reversal transformation 
T to ([9]) and using the anti-Hermiticity of T, we have 

Ulul = e- 2 ^c/t[/t. (11) 

The compatibility between (Q and (TlT!) leads to X a ^b = 0. 

In 3+1 dimensions, the only allowed statistics for particle excitations is boson or fermion, 
o-j = ±1. Then the permutation group representation is uniquely determined as rf = T a , 
with matrices T a satisfying 

T a T b = T b T a . (12) 

The commutation relations ([7]), ([H]) and Q now reduce to 

T a U b = e-^ i/3)5 ^U b T a , U a U b = e 27TiX ^U b U a , (13) 

with a, b = 1, 2, 3. The topological discrete algebra in the quark deconfmement phase consists 
of the flux insertion operators U a and the quark winding operators T a with the commutation 
relations (fl2l and (fl3"]). 

Now we show our main claim in this paper: // the ground state of QCD has a mass gap, 
the quark deconfinement phase has at least 3 3 -fold ground state degeneracy on T 3 . To show 
this, consider the following process. First, create a pair of quark and anti-quark out of a 
vacuum, then move the quark by T a . After the quark returns to the original position, we pair 
annihilate the quark and the anti-quark. Suppose that there is a mass gap to excitations 
above the vacuum space and these operations do not close the mass gap, then these processes 
define the operation of T a from a vacuum to a vacuum. Since T a (a = 1,2,3) commutes 
with each other, we can take the basis of the vacuum space which diagonalizes Ti, T2 and 
T 3 simultaneously, T a \r}) = e ir,a \r]),r) = (771, r] 2 , r] 3 ). By applying U a (a = 1,2,3) to this and 
using (fl3|) . we have 

T x {U[U°Ul\rj)) = e^-^^UlU^r)), 
T 2 (UlU^lrt)) = e lfe - 27rs/3) f/[£/ 2 s f^|r7), 

^3 {UlU^r))) = e^-^UZUlUllr,), (14) 
6 



where r, s and t are integers. Therefore, it is found that there are 3 3 distinct sets of 
eigenvalues of T a 's. This implies that the ground state (vacuum) in the quark deconfinement 
phase has at least 3 3 -fold degeneracy. 

On the other hand, if quarks are confined, the topological discrete algebra becomes trivial 
and the ground state degeneracy obtained above disappears: In the confinement phase, the 
permutation group for hadrons, not for quarks, classifies the physical states. The permuta- 
tion group for hadrons is also defined by ([6]) if <7j and t" are interpreted as those for hadrons. 
On the contrary to the quark deconfinement phase, however, all the generators of the per- 
mutation group for hadrons commute with the flux insertion operators U a (a = 1, 2.3). This 
is because hadrons have integer electric charges. From this property, the movement r" of a 
hadron around the inserted flux $ a = 2ir gives only the trivial Aharanov-Bohm phase, which 
leads to TfZ7& = U^Tf, (JiU a = U a <Ji,(a,b = 1,2,3). Then by the Schur's lemma, the flux 
insertion operator U a reduces to a unimodular constant for any irreducible representation 
of the permutation group for hadrons. In addition, since the allowed representations of the 
permutation group for hadrons are fermion and boson, rf for a hadron is again uniquely 
determined as rf = T a with mutually commuting matrices T a (a = 1,2,3). Because all 
the elements of the topological discrete algebra commute with each other, no ground state 
degeneracy is obtained in the confinement phase. 

Note that the ground state degeneracy obtained in the quark deconfinement phase de- 
pends on the topology of the space-manifold. This is easily seen by considering the system 
on a 3-dimensional box with the free boundary conditions, which is homotopic to the 3- 
dimensional ball B 3 . Because no non-contractable loop exists on this space-manifold, the 
permutation group consists of only the exchange operators <7j, and does not include r". 
Moreover, the operators U a (a = 1,2,3) can not be defined. So no topological discrete 
algebra is defined and no ground state degeneracy is obtained in this space-manifold. In 
general, if the space-manifold on which the system is defined has I independent spatial non- 
contractable loops, then the minimal ground state degeneracy in the deconfinement phase 
becomes 3 l . Our results here indicate that the deconfinement phase in QCD is topologically 
ordered, and the quark confinement and deconfinement transition is described properly by 
the concept of topological order. 

In the static limit of QCD, the minimal topological degeneracy obtained above is repro- 
duced by a conventional argument using the Wilson loop along the non-contractable loop 
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C a , W(C a ) = t^YlneCa U n ,a- In this limit, all the quarks are infinitely heavy and decoupled 
from the dynamics. So the system is effectively described by the pure 577(3) gauge theory 
The pure 577(3) gauge theory is invariant under the transformations, 

U n , a -> e 2mm/3 U n , a , (m = 1, 2, 3), n a fixed, (f 5) 

which rotate all spacelike links in the n a direction at a fixed n a by an element of the 
center of SU(3). The Wilson loop W(C a ) is transformed by this center symmetry as 
W{C a ) -> e i2mn / 3 W(C a ), so the expectation value (W{C a )) is an order parameter for the 
center symmetry. In the quark confinement phase, from the area law, it follows that in the 
temporal gauge (W(C a ,r)W^(C a ,7 J )) ~ e -t7jv °l r-T with the imaginary times r — it and 
t' = it'. Here a is a positive constant. Thus using the cluster property 

(W(C a ,r)W^C a ,r')) |T ^°° \(W{C a ))\\ (16) 

we have (W(C a )) = 0, so the center symmetry is not broken. On the other hand, in the 
quark deconfmement phase, it is possible that (W(C a )) ^ 0, and the center symmetry can 
be spontaneously broken. If (W(C a )) ^ for all C a 's (a = 1, 2, 3), we have 3 3 different set of 
(W(C a ))'8, which are related to each other by the center symmetry. Thus the ground state 
degeneracy is 3 3 -fold, and it coincides with the minimal ground state degeneracy obtained 
from the topological discrete algebra. 

The topological discrete algebra constructed above has a similarity to the 't Hooft algebra 

aa 

. For example our relation 

T a U b = e~ {2m/3)5 ^U b T a (17) 

corresponds to the following relation given by the 't Hooft, 

W{C)B{C) = B(C')W(C)e 27Tin/Af % M c = 3 (18) 

where C and C denote closed curves in 3-dimensional space, n the number of times the curve 
C winds around C in a certain direction, and B(C) the 't Hooft loop along C . However, 
there exist essential distinctions between them. First of all, the 't Hooft algebra can be 
defined only when the dynamical quarks are absent, but our topological discrete algebra is 
defined in the presence of the dynamical quarks. Second, the topological discrete algebra in 
the quark deconfmement phase is different from that in the quark confinement phase, but 

8 



the 't Hooft algebra is the same in any phases. Third, new quantum numbers X ajb , which 
are missing in the 't Hooft algebra, exist in the topological discrete algebra. 

Now we would like to address the physical meaning of X a>b . Consider the degenerate 
ground states <pK — \$)k {K = 1, • • ■ , d) in the presence of inserted fluxes <fr = ($1, $2, ^3)- 
They satisfy U a \<f?)x — e* 7a< -*' ) |$ + d2ii)x, where 7a ( < £a) is the quantum holonomy given by 

From U a U b = e 2niXa > b U b U a , it follows 

7a ($ + 62tt) + 76 (#) - 76 ($ + a2vr) - 7a (*) = 27rA a , 6 + 2vrM, (20) 

where M is an integer, which is independent of K because in general the degenerate ground 
states are related to each other by some symmetry. Then by the Stokes's theorem, the Hall 
conductance a ab (a ^ b) [2^] 



e 2 d f 2n f 2n 



hd ^ Jo Jo 2ni 



(21) 



reduces to 



e 2 



a ab = j{X a ,b + M). (22) 

This indicates clearly that a fractional \ a , b implies the fractional quantum Hall effect. 

The above discussions can be generalized to the SU(J\f c ) QCD. For this purpose, it is 
convenient to introduce a fictitious U(l) gauge field as an external field, instead of the 
electromagnetic gauge field, and assign the fictitious U(l) charge Q = 1/M C to all quarks. 
Then, consider an adiabatic insertion of the fictitious U(l) flux by 2tt through the hole 
h a of T 3 , which is represented by a unitary operator U a (a = 1,2,3). The quark winding 
operator T a is defined in the same way as above. The topological discrete algebra in the 
quark deconfinement phase is given by 

T a T b = T b T a , T a U b = eW^»W a , U a U b = e 2 ^ /l ^UyU a , (23) 

where k a ^ and l a)b are co-prime integers, and l a ,b is a divisor of Ac- Using these relations, 
the minimal ground state degeneracy on T 3 in the quark deconfinement phase is found to 
be A/" 3 . On the other hand, in the quark confinement phase, the topological discrete algebra 



is trivial and no topological degeneracy arises because the charge Q of any hadron is an 
integer. 

To count the ground state degeneracy, we have assumed that the system has a finite gap. 
While color charges are screened in the presence of the gluon mass gap, our result indicates 
that the quark confinement is not synonymous with the color screening. In addition, the 
topological discrete algebra itself is valid even in the gapless system. The concept of ground 
state degeneracy becomes subtle in the gapless system, but the degeneracy could be identified 
by examining the finite scaling carefully. 

To conclude, we have argued phases in QCD by using a discrete symmetry algebra which 
is manifest only on a space with non-trivial topology. The topological degeneracy of the 
ground state, which indicate the presence of a topological order, is derived and it is found 
that even in the presence of dynamical quarks it is a good quantum number distinguishing 
the quark confinement phase from the deconfinement one. 

The author would like to thank H. Kawamura for helpful discussions. 
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